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In this paper, we present a model for the effects of the tachochne on the dif- 
ferential rotation in the solar convection zone. The mathematical technique relies 
CN '. on the assumption that entropy is nearly constant ( "well-mixed" ) in isorotation 

surfaces both outside and within the tachochne. The resulting solutions exhibit 
Pm ! nontrivial features that strikingly resemble the true tachochne isorotation con- 

tours in unexpected detail. This strengthens the mathematical premises of the 
P_i! theory. The observed rotation pattern in the tachochne shows strong quadupolar 

Q ! structure, an important feature that is explicitly used in constructing our solu- 

tions. The tachochne is treated locally as an interior boundary layer of small but 
ci '. finite thickness, and an explicit global solution is then constructed. A dynamical 

link can thus be established between the internal jump in the angular velocity 
►^ , at the tachochne and the spread of angular velocities observed near the solar 

O ! surface. In general, our results suggest that the bulk of the solar convection zone 

is in thermal wind balance, and that simple quadrupolar stresses, local in radius, 
•^ I mediate the tachochne transition from differential rotation to uniform rotation 

\^ ' in the radiative interior. 

o ■ 
o 

Subject headings: convection — hydrodynamics — stars: rotation — Sun: rota- 
tion — Sun: helioseismology 
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1. Introduction 

Recent work suggests that the isorotation contours in the bulk of the solar convection 
zone (hereafter SCZ), correspond to the characteristic curves of the vorticity equation in 
its "thermal wind" form (Balbus 2009, hereafter B09; Balbus et al. 2009, hereafter BBLW; 
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Balbus & Weiss 2010). In its undeveloped form, this equation relates entropy gradients 
to angular velocity gradients. The mathematical existence of the characteristics demands, 
therefore, some sort of functional relationship between the entropy S and angular velocity 
Q. The theory presented in the above papers is predicated upon the notion that within 
constant fl surfaces, the dominant, long-lived convective cells carry out their task of mixing 
entropy with great efficiency. This means that within a constant Q surface, the excess 
entropy gradient above and beyond the minimal radial gradient needed to maintain the 
convective state is, in essence, eliminated. The remaining "residual" entropy is thus uniform 
in constant Q surfaces. More explicitly, the residual entropy S' = S — Sr is formed by 
subtracting a function of spherical radius r (something close to the driving background 
entropy that has been averaged over angles), and S' is set equal to a function of Q. This is 
a powerful constraint, providing the necessary mathematical link for interpreting the solar 
isorotation contours as trajectory characteristics of a quasilinear partial differential equation. 
BBLW advances several different lines of argument (heuristic, numeric, and goodness-of-fit) 
supportive of this picture. Originally an MHD basis for the S — fi connection was put forth 
(B09), but the greater simplicity and generality, as well as the broad agreement with many 
different lines of argument, together suggest primacy for the BBLW mechanism. 

An attractive feature of the theory presented in BBLW is that it provides simple and 
intuitive answers to not-so-trivial questions. For example, why can't the solar convection 
zone rotate on cylinders? Because if it did, then Q would necessarily depend upon the 
spherical colatitude angle 6, therefore the entropy would also depend upon 6, and this is 
inconsistent with a vanishing z gradient in the thermal wind equation (TWE; cf. equation 
[2] below). Why then is Q dominated by its 6 gradient? Now we argue the other way: 
because the residual entropy has its dominant radial gradient removed (physically by mixing, 
mathematically by a simple subtraction), the relative importance of its 6 gradient has been 
enhanced. But residual entropy shares the same isosurfaces as Q. In numerical simulations, 
the fact that imposed latitudinal entropy gradients enforced at an interior boundary lead 
also to latitudinal gradients in Q (Miesch, Brun, & Toomre 2006) becomes tautological: 
latitudinal S'-gradients lead inevitably to latitudinal fi-gradients because these gradients are 
always counteraligned. 

In its current formulation, BBLW theory addresses the differential rotation in the Sun 
away from the tachocline. This is nothing if not prudent: the dynamics of the solar tachocline 
is a notoriously controversial and difficult topic (e.g. Hughes, Rosner, & Weiss 2007). Yet, 
in studying the isorotation contours within the tachocline, one is struck both by the highly 
localized nature of the distortion, and by the relative simplicity and uniformity of the their 
appearance. In contrast to the bulk of the convection zone and the outer boundary layer, 
there is a sense of inevitability to these curves. Perhaps their gross form might not be 
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sensitive to the complex details of the internal dynamical structure of the tachocline, any 
more than the gross adiabatic temperature profile of the SCZ depends upon the details of 
the convective heat flux. We might look instead to the exigencies of local forcing and to 
the strikingly quadrupolar structure of the region. The resulting mathematical demands 
of joining the isorotation contours in such a forced region to the overlying convective zone 
contours treated in BBLW may be sufficiently restrictive that the role of turbulence is tightly 
confined. In its current form, BBLW theory has too much "free play:" it needs to be 
restrained. In this work, we will develop this point of view, and pursue its consequences. 



2. Analysis 

2.1. Trajectory and solution characteristics of the TWE 

We follow the notation convention of BBLW. Cylindrical coordinates are denoted by 
radius R, azimuthal angle and vertical coordinate z. The spherical coordinates are (r, 6, 0), 
where r is the radius from the origin, 6 the colatitude angle, and the azimuthal angle. The 
angular velocity Q, pressure P and density p are understood to be azimuthal averages. The 
dimensionless entropy function a is defined by: 

a = \nPp-"', (1) 

where 7 is the usual specific heat ratio, or adiabatic index. 

In SCZ applications, the TWE is the dominant balance of the vorticity equation, after 
azimuthal averaging, between the large scale rotation and the baroclinic terms. The latter 
arises from a mismatch between equipotential (or isobaric) and isochoric surfaces. More 
specifically, contributions from convective turbulence are ignored (i.e., the convective Rossby 
number is small [Miesch & Toomre 2009]) as are those from from magnetic fields. The TWE 
for the convective zone is (e.g. Thompson et al. 2003, Miesch 2005, B09): 

n'-f = ^|, (2) 

oz 7r 06 

where g = GMq/t'^, G is the gravitational constant and Mq is a solar mass. In BBLW 
theory, there is a functional relationship of the form 

a' = a-ar = /(fi^) (3) 

where a^ is any function of r (in practice something very close to an angle-averaged a) and 
/ is an unspecified function. Replacing a' by a does not alter the TWE, and with the aid of 
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equation ([3]), the TWE ([2]) becomes 



dVi^ ( gf tan^\ dQ^ ^ 

' + Ktt = 0> (4) 



dr \7r^ sin 6' cos 6* r J 09 

where /' = df/dQ'^. In this study, we wish to retain, at least in a formal sense, the tachocline 
meridional stress T{r, 6) on the right side of the equation B 

^^' ( 9f . tane\ on' 

T(r,6). (5) 



(7) 



dr \7r2 sin 6* cos ^ r J 86 

Although T is certainly poorly understood in detail, there is presumably some function of 
position that makes equation ([5]) correct. The point is that the ensuing formal solution is 
not without content. The angular velocity Q^ satisfies 

^^T(.,«, (6) 

on 

d9 _ _ gf'in^) _tan9 

dr 7r^ sin ^ cos r 

In other words, to solve for Q{r,6), we use the contours as defined by the characteristic 
differential equation of BBLW theory; along these contours Q changes according to equation 



We shall assume that equation ([3]) is valid both in the bulk of the convective zone as well 
as in the tachocline. This is by no means obvious. The low Rossby number approximation 
is likely to hold well in the tachocline and convective rolls should be strongly affected by 
shear. But since the turbulent stresses in the tachocline are very different from those in the 
convective zone, even if such a functional relationship holds, the function / may in principle 
be quite different in the two regions. 

On the other hand everything changes in the two regions, not the least of which is 
the sudden appearance of the stresses themselves. For the class of solutions we consider 
here (which reproduce the data strikingly well), a change in /' cannot be disentangled from 
a change in T, a point we will explicitly discuss (§2.6). Hence, we retain the functional 
equation ([3]) throughout the convective zone and the tachocline. 



^ In passing from the lower convective zone into the radiative layers, the baroclinic term proportional 
to {dP/d9){da/dr) should be formally be retained in our equation. As a practical matter, the rotational 
pattern is essentially spherical in this region, and the dominant balance is between the inertial terms and 
the external stress. (We thank the referee M. Miesch for drawing our attention to this point.) 
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Note that, in contrast to the vanishing T hmit, when T is finite /'(f2^) is generally 
not constant along the characteristic. But this does not mean that there is nothing to be 
learned from considering the constant ^2 as a special case for the current, more general, 
problem. Indeed, in our first approach we should keep formalism to a minimum. It has 
been earher noted (B09) that choosing fi^t) to be a global constant captures important 
semi-quantitative features of the true SCZ isorotation contours, so let us begin first with 
this mathematically simple case. We will then be poised to look at a more complex problem. 



2.2. /'(fi2) = constant 

In the limit of constant /', the formal trajectory characteristics are identical to those in 
BBLW theory, both outside and within the tachocline: 



sm 9 = ^ sm fc^o H ^ 



where tq is the radius at which Q is initially specified (generally near or at the surface), and 
^0 marks the beginning value of 6 for a particular characteristic path. The sign of T is now 
clear: since along a BBLW characteristic, Q increases with decreasing r, T must be negative 
at high latitudes and positive at equatorial latitudes. Along such a characteristic, denoted 
6{r), Q'^ is given by 

pr 

n^ = nl+ / T[r,e{r)]dr (9) 

Jro 

where Qq is the value of Q at the defining surface. (Note that ro, which is generally taken 
at or near the surface, may be larger than r.) The formal solution (Q is in fact general for 
any /', provided that 6{r) is taken to be the proper trajectory characteristic. 

It is not difficult to extract the qualitative behavior of Q from equations ([8]) and iQ. The 
stress T is very small except in a very narrow region near the transition radius tt- Thus fl^ 
remains essentially fixed at its initial value ^q along most of the extent of the characteristic. 
This, of course, is BBLW theory. Then, there is a sudden increase in T very near r^, and 
along the BBLW characteristic 6{r), Q now makes a very sharp rise (or a sharp drop near 
the equator) as it settles to the uniform rotation rate in the radiative interior. 

The resulting isorotational curves can be understood with the aid of figure (1). This 
diagram is a local representation of the SCZ-tachocline boundary in a region where the 
angular velocity begins to increase sharply moving inward. The radius r increases upward 
to the right; 6 increases downward to the right. The numbers 1' though 4' label BBLW 
characteristics, which in the convection zone are precisely the isorotation contours. Hence 
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these numbers may also be thought of as angular velocities, increasing with the magnitude 
of the number. 

Consider the rightmost contour labeled 4'. Between 4' and the dot labeled 4 on the same 
curve, the contour passes through the convective zone, and the angular velocity remains 
constant. Dot 4 marks the beginning of the tachocline, and a steep rise in the angular 
velocity. To continue the trace of the isorotation contour, a sharp turn to the (reader's) left 
is required, and the curve continues by passing through dot 4 on the curves labeled 3', 2', 
and 1', etc. Thus, within the narrow radial band of changing Q, if one were to follow curves 
of constant Q, from one BBLW trajectory characteristic to those immediately adjacent, a 
quasi-spherical (r increases somewhat with 6) path would emerge. This abrupt change in 
the isorotation contours is clearly seen in the helioseismology data. The more delicate 6 
dependence of r is also present in the data, though this may be an artifact of the inversion 
technique. 



2.3. Green's Function Solution 

The above description motivates a Green's function approach. Consider a tachocline 
stress of the form 

T{r, 6) = -{An^)F{e)5{r - r^) (10) 

where AQ^ is the "jump amplitude," F{6) is a function of 6 that modulates the magnitude 
of the radial inward jump (F is positive at high latitudes, negative at small latitudes), and 
6{r — rj) is the Dirac delta function. Then along a BBLW characteristic, 

n'' = nl + Q{rT - r){An^)F{e{rT)) (11) 

where 9{r) may be determined explicitly from equation (|8]), and G(a;) is the Heaviside func- 
tion. (Recall that 6 is unity when its argument is positive and zero when its argument is 
negative.) In this rather drastic idealization, BBLW theory breaks down at a single radius 
tti which we identify with the tachocline. The BBLW isorotation contours change abruptly 
to a single spherical shell in this model, compressing the entire tachocline down to a range 
of vanishing thickness. 



2.4. Quadrupole Structure 

It is possible, of course, to continue along the lines of the previous section by introducing 
progressively more sophisticated models for T{r,6) and f'{Vl?)^ and in fact we shall shortly 
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consider the case in which /' is a hnear function of Q^. It is aheady clear, however, that 
within this scheme, the hehoseismology data can be reproduced — at least qualitatively. If 
our approach is not quite tight enough to be broadly predictive, it might be more profitably 
applied in the reverse sense: use the solar rotation data directly to infer the form of T{r, 9). 

A very important result emerges immediately. In figure (2), on the left we see the isoro- 
tation contours taken from GONG dataci, together with the BBLW characteristics overlaid 
for reference. As one passes through the tachocline from larger to smaller radii, the inward 
jump in Q"^ changes from positive to negative for 6 larger than 54.7°, the value of 6 for which 
cos^ 6' = 1/3. The inset to the right in figure (2), displaying GONG data together with the 
location of the critical angle, shows this very clearly. The magnitude of the jump then grows 
(in the negative sense) as 6 increases, as we move toward the equator. The spacing of the 
isorotational contours is, at least crudely, about half the interval at the poles as compared 
with the equator. Even granting the inevitable uncertainties associated with the polar data 
(which are not reliable poleward of 70 degrees), these results are much what one would expect 
if F{6), the angular component in our radially local model, were proportional to 

P2(cos^) = -(Scos^^-l), 

the Legendre polynomial of order 2. To leading order, the tachocline stresses are revealing a 
dominant quadrupolar structure. 

This is such a fundamental property of tachocline structure, it is curious that the near 
coincidence of the zero of P2(cos 9) and the "zero" of the tachocline has not been emphasized 
in the earlier literature. While meridional flows and their associated stresses have been 
an important feature of dynamical models of the tachocline (e.g. Spiegel & Zahn 1992, 
Gough & Mclntyre 1998), and such flows are generally associated with P2{cos9) stresses 
(Schwarzschild 1958), what is striking here is the manifest quadrupolar morphology. In §2.6 
below we use explicit quadrupolar forcing to leverage a mathematically precise isorotational 
contour solution from our fundamental equation, with essentially no other input. So it is of 
some interest to try to understand its origin. 

The most important omission from the thermal wind equation, the basis of our analysis, 
is the contribution from turbulent stresses and meridional flow. These may be kinetic or 
magnetic, but they will generally take the form of a sum of terms proportional to 

di{ViWjR) 



^We thank R. Howe for providing us with these results. 





Fig. 1. — Schematic local rendering of the isorotation contours in the presence of tachocline 
meridional stress. CZ denotes the convective zone and T denotes the tachocline. The thick 
curves labeled 1', 2', etc. are the CZ isorotation contours in BBLW theory. Numbers may 
be thought of as representing the relative ordering of fiducial Q values. The presence of the 
tachocline, indicated by black dots, is to increase fl along what was an isorotation contour. 
The new, proper isorotation contours (thin lines) resemble those seen in the helioseismolgy 
data: quasi-spherical, with an increase in r as 6* increases, and with a sharply angled upward 
turn toward the surface. 





P.(cosfi) = 0, 8 = 54.7° 



Fig. 2. — Left: Constant Q contours (white) of the thermal wind equation from BBLW 
theory overlaid on top of (black) isorotation contours (GONG data, courtesy of R. Howe). 
Blue curve is the bottom edge of the convective zone. Right: Tachocline is absent when 
P2{cos9) = 0, consistent with quadrupole structure. 
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where the index i is itself summed over, and V and W represent symbolic vector quantities 
(velocity fluctuations, vorticity, current density, etc.). It is the r and 6 components of i 
that are important (the (p components vanish either explicitly or upon azimuthal averaging); 
hence the earlier designation of the resulting stress as "meridional." The net contribution of 
these correlated products evidently behaves as a large scale quadrupolar forcing. Given the 
symmetry of our flow, it is not surprising that the leading order departure from spherical 
structure is quadrupolar, but it is rather surprising how well this approximation works. 



2.5. F(9) = P2{cos9): an explicit solution 

We now combine the results of the previous section to obtain an explicit mathematical 
solution for Q. Equation ([8]) implies quite generally that 

1 1 \ 

fl2) 





sm 6^0 = ^ sm 6 ^ | | . 




rl 


rol \r roj 


Now at r : 


= tt equation ^ gives. 






sin^ 9{rT) - "'" 


2 , 2/'GMo / 1 1 
sm^9o+ ■' , Q 



(13) 



Then, substituting equation flT2|) in flT3l) leads to 



sm 9{rT) = — sm 9 -\ ^ ( J , (14) 



2 n/ \ ' -2 

rj. r^7 \rT r 

or 

cos 9{rT) =1 ^ sm 9 ^ ^ , (15) 

\ r|,/ r|,7 \r rj-J 

Then, equation flTTl) becomes 

fi2(r,e) =f]g(sin2 0o) + 0(rT-r)(Afi2)P2[cos^(rr)] (16) 

where sin^ ^q is given by equation flT2l) . cos 6'(rT) by equation (IT5|) . and we have assumed 
that F is the Legendre polynomial of order 2. 

Equation ( TT6l) is the desired explicit solution for i7^, but there is another way to formulate 
the content of our simple problem: given Ail^, what kind of a spread of surface velocities is 
associated with this jump? The answer is 

fi2(sin2 ^o) = ^L - (Afi2)P2[cose(rr)] = fi^ _ (^^2) ^ 3^^^2) sin2[^(^^)] (17) 
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where flrad is the angular velocity of the uniformly rotating radiative interior. With the aid 
of equation flT3l) . we may rewrite this in terms of sin^ 9q. Choosing VQ/r"^ = 2, and denoting 
the final term of equation flT^ as C, (a number of order, but typically less than, unity), we 
find 

nlism-" Oo) = fiL - (1 + 0(A^') + 3(Afi2) sin2(^o) (18) 

This calculation is a bit too simple to generate more complexity than a sin^ Oq departure 
from uniform surface rotation. Nevertheless, for reasonable values of AQ^, equation (fT8l) 
gives a very respectable rendering of a solar-like surface rotation profile. More importantly, 
it establishes a dynamical coupling between the spread in Q^ present near the surface and 
the jump in Q^ at the location of tachocline. The connection between these two regions is 
possible because of the existence of linking trajectory characteristics. 

Equation (ITSl) should not be construed as a statement that the spread in surface angular 
velocities is actually caused by forcing from the tachocline; indeed, the sense of causality is 
generally thought to run in the opposite sense: from the SCZ to the tachocline (e.g. Gough 
2010). Instead, equation (fTSjl shows explicitly how the dynamics of thermal wind balance 
plus external driving leads to a relatively simply coupling between the radial jump in Q^ at 
the tachocline and the angular spread in Q^ at the surface. 



2.6. f'{Q^) linear in fi^ 

2.6.1. Governing equation 

We next consider the case in which /' varies linearly with Vl? . This is useful both in 
establishing that there is no particular sensitivity to the assumption of a constant /', and in 
allowing more detailed structure to be revealed. In fact, the essential qualitative effects of a 
nonconstant /' will be evident once the governing equation is at hand. 

We begin by rewriting the trajectory characteristic equation (jTj) in the compact form 

where, as before, R = r sin 6, the axial radius. Let 

f'm = -« - (3n', (20) 

where a and (3 are constants. Dividing by g, differentiating with respect to r, and using 
equations (1201) and ([6]) leads to 

y-'-f-'-^nr,0> (21) 

ar g ar 7 
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Recalling that g = GMq/v'^, this leads to the final form of our equation for the trajectory 
characteristics, now free of all explicit Q dependence: 

We now make the very reasonable and important assumption that T{r, 9) is local in r, so 
that this variable is always near the shell r = tt whenever T is finite. In essence, we ignore 
curvature terms, something we must do self-consistently since these global terms are not 
known for general T{r, 9). The local model avoids this difficulty while retaining the essential 
physics, because the tachocline transition is not, in reality, very broad. 

Consider equation ( l22l) when the driving term T{r, 9) takes the (radially local) form 

T(r,0) = -ToP2(cos0) (23) 

where Tq is a constant term. We assume that T is present only for r < tt- In this region, 
the differential equation f l22|) becomes 



—5 = - sm t^ — 1 



(24) 



It should be noted that the function /' makes its appearance here only through the (3 
parameter, and then only via the combination PTq. This illustrates the point first made in 
the Introduction: a different value of /3 can equivalently be regarded as keeping /3 the same 
and changing Tq. There is nothing to be gained by allowing a change in /' in the tachocline, 
at least not for the class of model of interest here. 

Replacing r by r^ everywhere leads to the differential equation 



dr'^ 7r|, \ 3/ \ 3 

This is a linear second order equation for sin^ 9, and the solution is 

2 
sin^ 9 = Ai cosh(A;r) + A2 sinh(A;r) H — (26) 

o 

where Ai and A2 are integration constants to be determineqj- 



•^ Equation (1241) may in fact be solved analytically without making the local approximation. But since 
locality has been used implicitly in ignoring the radial structure in T, it is best to be self-consistent. 
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2.6.2. Global solution 
The characteristic solution for r > r^ is aheady known (BBLW): 

W« = 4w«„ + ?5^^^fi-i) (27) 

where 

ng(sin2 ^o) = fi? + f^2 sin^ ^o, (28) 

and Ql and Ql are fixed constants. 

We seek global solutions in which the exterior solution (!27|) for r > tt joins smoothly to 
the interior solution (126 p for r < rj' at r = r^. The solution for sin^ 9 and its first derivative 
with respect to r should be continuous. Two dimensionless parameters are needed, defined 
by 

7r;f 

^ 2GM,(« + ^n;)^ 

An easy way to determine representative values of .^i and .^2 is to use the B parameter 
introduced in B09 for solving for the form of the SCZ isorotation curves: 

_B^ 2GMrj'{nl)nl ^g^^ 

(Note that the solar radius notation r© was used in B09 instead of our tq.) A very simple 
and reasonable SCZ model is discussed in B09 using 

-^ = 0.12 + 0.8 sin^^o, (32) 



which would imply 



ft=(Ii)'o.8, 6=(^)'o.l2. (33) 



Two other parameter combinations may be defined that will be later useful: 



6 = 1^ 



^4 = 6 



2 



^0 






6 (34) 

(35) 
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At r = tt, the second derivative of sin^ 9 along a characteristic is discontinuous, but its 
value and first derivative are continuous. With equations fl26|) and fj27|) . together with the ^ 
parameters, these two continuity conditions may be written: 

2 

,^3 sin^ ^0 ~ ^4 = ^1 cosh krT + A2 sinh krx H — , (36) 

(^1 - 2^3) sin^ 6^0 + 6 + 2^4 = krriAi sinh krr + A2 cosh krr) (37) 

Solving for Ai and A2, we find the following interior solution for the characteristics: 

sin^g-^=[e3sin^go-e4-2/3]coshA:(rT-r)- '^^^^^^^^J~''^ [(6 - 2^3) sin^ go + 6 + 2^4] 

(38) 
For future reference, we list here the form of this characteristic isolating sin^ ^o^ 

. 2n sinH- 2/3 + {^i + 2/3) cosh k{rT-r) + sinh k{rT-r){^2 + 2^i)/krT ,^^. 

sm uq = — —, r : — - — —, . , . . (39) 

^3 cosh k[rT — r) + smh k{rT — r){2t,3 — q,i)/ kr^ 

To solve for Vt^ along the characteristic ( l38l) . use equation (jQ]), noting that T{r,9) van- 
ishes for rQ > r > vt'- 

T{r, 9)dr = Ql + ^^ (sin^ 9{r) - 2/3)dr (40) 



Integrating with the help of equation ( l38l) we obtain 

3T 

fi2 = ^2 ^ —^ [((^, _ 2^3) sin^ ^0 + 6 + 2^4) (l/A;VT)(cosh k{rT - r) - 1) 

-(^3 sin^ ^0 - 6 - 2/3)(l/fc) sinh k{rT - r)] (41) 

This gives us the solution along the characteristic beginning at the surface value r2g, which is 
itself a function of sin^ ^o- We would like to have an expression for Vt^{r, 9), without reference 
to the characteristics, so that we may obtain the isorotation contours. To this end, we first 
use equation ( 128|1 to expand fi^, followed by equation (139|) to eliminate sin^ 6*0 from the above 
expression. In this way, we may obtain the desired isorotation contours in the form of sin^ 9 
as a function of Q^ and r (as well the various parameters of course). The result of this rather 
lengthy algebraic excursion is that the isorotation contours in the tachocline are given by: 

• 2/) 2 Hi - H2 + Hs 



where 

fi^ — fi? r^ , ,, , , ^ ^ , sinhA;(rr — r) 

^ ^,coshk{rT-r) + (2^3-^1)- ^ ^ > 



Hi 



Ql 



kri 



(43) 
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H2 = [(^4 + 2/3) cosh k{rT - r) + (^2 + 2^^){1 /krr) sinh k{rR - r)] (44) 

Hs = [6(6 + 26) - (6 + 2/3)(26 - 6)] """'^^^"^'"^"^ (45) 

rr kvT (, . ^,, .^ (26-6)(coshfc(rT-r)-l) ^^ 

Ha = -y- I 6 smh k{rR - r) ^ I (46) 

Equations (138|) and f l42|) for the characteristics and isorotation contours respectively comprise 
the interior tachochne solution for Q{r,6). 

Notice that we have not made use of any inner boundary condition enforcing uniform 
rotation at a particular location marking the beginning of the radiative zone. In principle, 
this location is part of what determines the form of the T{r, 9) stress. In practice, however, 
the simple model that we have adopted reproduces the helioseismology data so well that we 
may we just stop the calculation at a value of r near 0.7. The isorotation contours at this 
point are very nearly spherical (cf. figure [4]). 



2.6.3. Modeling the Sun 

At a transition radius of tt = 0.77ro, the helioseismology data show a modification of 
what we refer to in this work as the "exterior" solution. This translates to 

6 = 1.75233, 6 = 0.26285, 6 = 1-28539, 6 = 0.060456 (47) 

The quantity {Q"^ — fiD/fil is simply sin^ 6^0, where 9o is the surface colatitude angle of the 
isorotation curve. Plots of the resulting characteristics and isorotation contours are shown 
in figures (3) and (4) for the case krx = 3. 

These figures are remarkable. The characteristics in figure (3) are everywhere smooth 
and regular. They are identical to isorotation contours for r > r-r, and bear no resemblance 
to them for r < r^- The corresponding isorotation contours (figure [4]) are continuous, but 
have sharply discontinuous first derivatives at r = r-p. What is astonishing, however, is 
the extent to which they resemble the helioseismology GONG data (cf. figure [2]). Recall 
that our model consists, in its entirety, of i) the assumption of cospatial surfaces of constant 
residual entropy and angular velocity in the thermal wind equation; and ii) local spatial 
forcing via an inhomogeneous term proportional to P2{cos6) that turns on at r = r^-. There 
is nothing else. Indeed, to keep the mathematics as simple as possible, we have not even used 
the best fit for the exterior solution shown in figure (2), but rather something much more 
basic. Still, the similarity of our closed form solution to the GONG data is unmistakable. 
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Fig. 3. — Trajectory characteristics (eq. [38]) for the SCZ and tachochne. The transition 
radius, shown as a dotted hne, is 0.77 of the surface radius, and we have chosen fcrj- = 3. The 
curves, which are here displayed well beyond their relevant physical domain, are everywhere 
regular and have continuous first derivatives at r = ry. 
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Fig. 4. — Isorotation contours (eq. 112]), with r^ = 0.77ro, krx = 3. The curves are 
continuous at r^ but have discontinuous first derivatives. The ensemble of isorotational 
contours bears a clear and striking resemblance to the helioseismology GONG data (cf. 
figure [2]). 
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Our results can be compactly summarized: the analysis shows that the SCZ is essentially 
in thermal wind balance from near the surface down to a transition radius of tt = 0.77ro, 
at which point distinctly quadrupolar tacholine meridional stresses bring the rotation profile 
to that of a solid body. For all of its detailed complexity, the role of turbulent angular 
momentum transport in the SCZ is strictly constrained, much as thermal transport has 
little choice but to establish a nearly adiabatic temperature profile. 



3. Conclusion 

The level of dynamical detail that can now be elucidated in numerical simulations of 
the SCZ is truly impressive (e.g. Miesch et al. 2008), but simulations incorporating the 
tachocline into global models of solar differential rotation are only just beginning to take 
form (Browning et al. 2006, Miesch et al. 2009). The solar tachocline remains a complex, 
poorly understood region of turbulent activity. While the details of the interaction between 
the tachocline and the large scale differential rotation of the SCZ are not yet within reach, 
certain gross features of this interaction may in fact be tractable. In this work, we have made 
an inroad by assuming that entropy is well-mixed in surfaces of constant angular velocity 
within the convectively unstable portion of the tachocline, just as it appears to be within 
the SCZ. Since the convective Rossby number is likely to be even smaller here than in the 
SCZ, this would seem to be a reasonable approximation. 

The explicit separation of the more general tachocline stresses from the mean inertial 
and baroclinic terms in the governing equation of vorticity conservation in its developed 
thermal wind form, and their reduction to a single formal inhomogeneous forcing term, has 
proven to be mathematically advantageous. The precise nature of this stress is, of course, of 
an unknown character. We would argue, however, that its consequences are not unknowable. 
The stress operates locally in radius and exhibits a definite quadrupolar angular structure: 
this is a very important point we wish to emphasize. It can be put to good use, because the 
formalism of quasilinear partial differential equations may be brought to bear on the problem. 
In particular, one may analyze the problem in terms of trajectory characteristics, whose form 
depends only implicitly on the tachocline stress terms, and solution characteristics, along 
which Q changes explicitly due to these same terms. 

In BBLW theory, the solution characteristics convey constant Q along the trajectory 
characteristics. These characteristics are then identical to the isorotation contours. It is 
thus of considerable interest that in a simple and revealing class of solutions, the trajec- 
tory characteristics are mathematically identical to those of BBLW theory, even within the 
tachocline. The solution characteristics are certainly altered, however: they dictate that Q 
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is initially constant along the trajectories, but then changes rapidly along these same curves 
as the tachocline is penetrated from above. The result is quasi-spherical isorotation contours 
inside of the tachocline, followed by a sharp upturn when the convection zone is reached from 
below. This behavior is also clearly seen in a more sophisticated, but still explicit, solution 
that is presented in this paper (cf. figure [4]). 

The radially local nature of the tachocline stresses is the crux of the problem, as it al- 
lows for mathematical modeling independent of the need of detailed knowledge of the spatial 
structure of the forcing. The simplest approach is to use Green function techniques with 
forcing by a Dirac (5- function at a transition radius tt- The dominant P2(cos^) angular 
structure of the tachocline stresses becomes very apparent when the radial behavior is sim- 
plified in this way. Since the monopole response is evidently small, a reasonable surmise is 
that the turbulent stresses in the tachocline are driven by off-diagonal terms in a correlation 
tensor. 

Angular momentum stresses within the SCZ create differential rotation, even if the 
initial condition is one of uniform rotation. Thus, in a causal sense, the tachocline likely 
owes its existence to the stresses of the convective zone, and not vice-versa. In its current 
state of development, our calculation cannot directly address this causal link, but it does 
suggest how a dynamical couple arises between the surface spread in angular velocity and the 
tachocline angular velocity jump. This characteristic-based connection between the jump in 
the angular velocity at the tachocline and the spread of angular velocities at the solar surface 
is an important consequence of simple quadrupolar forcing. 

Venturing even a little beyond a simple Green's function approach gives a much richer 
return. If one allows both an f2-dependence in our /' function and a finite thickness for the 
transition region, the mathematical structure of the problem becomes more complicated: the 
solution and trajectory characteristics are now coupled. But it is still possible to obtain an 
explicit closed form solution for the isorotation contours. This, the principal result of this 
paper, is shown in figure (4). The relatively simple solution so much resembles the GONG 
contours seen in figure (2), there can be little doubt that at the very least, the mathematical 
content of our approach is sound. Even if viewed at a purely phenomenological level, this is 
surely progress. 

While our results provide a simple framework to aid in the elucidation of the global 
character of solar differential rotation, many questions remain unanswered. Quite apart 
from the grand challenge problem of tachocline structure and dynamics, it is legitimate to 
question the work presented here on its own terms. How general and robust, for example, is 
the separation embodied in equation (5), or the critical assumption that a functional relation 
remains between residual entropy and angular momentum in the lower buoyancy tachocline? 
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Why should such a simple approach work so well? There is also evidence that the tachocline 
thickness varies with 9 (Christensen-Dalsgaard & Thompson 2007), an effect completely 
neglected here. The relative simplicity of our arguments is somewhat deceptive, because it 
masks the deeper conceptual problem of understanding precisely how the coupling between 
the tachocline and the convective zone rotation profile is established (see, e.g. Rempel 2005, 
Gough 2010). These issues certainly have not been addressed by this work; doing so without 
recourse to turbulence modeling will remain the province of numerical simulation for the 
foreseeable future. 

But in this weakness also lies some power. We have noted earlier that the SCZ adiabatic 
temperature profile sheds little light on the dynamics of convective transport, though it is in 
an important sense completely beholden to it. This is not usually viewed as a shortcoming; 
rather it is a useful and powerful constraint. This is a point worth remembering: some 
aspects of the solar rotation problem are bound to be less mysterious than others. If thermal 
wind balance and the confluence of isorotation and constant residual entropy surfaces lead 
to important constraints for the Sun's rotation profile, the theorist's task will be greatly 
simplified — and understanding how the Sun rotates may prove to be a less daunting challenge 
than it might have first appeared. 
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